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ABSTRACT 


A  digital  data  acquisition  system  was  implemented  for  the  study  and  collection 
of  scattered  waveforms  from  various  simple  objects.  Reflected  waveforms  were 
extracted  from  interference  with  other  signals  by  a  simple  subtraction  technique. 
Backscattered  and  bistatic  reflections  were  observed  and  recorded.  Various  effects 
on  scattered  waveforms,  including  target  composition,  pulse  length,  and  frequency, 
were  experimentally  observed.  Bistatic  scattering  results  from  a  stryofoam  sphere 
at  a  ka  value  of  10  were  in  agreement  with  theoretical  predictions  for  a  rigid  sphere. 
Bistadcally  measured  data  at  70  kHz  from  a  styrofoam  cylinder  of  length  20  cm  and 
diameter  2.5  cm  followed  a  trend  for  bistatic  angles  up  to  150  degrees  but  became 
amorphous  as  the  bistatic  angle  approached  180  degrees. 
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Any  object  with  a  characteristic  impedance  different  from  the  surrounding 
fluid  will  affect  an  acoustic  wave  impinging  on  it.  The  wave  will  be  transmitted 
through  the  object  and  scattered  from  it.  Flexural  waves,  compressional  waves,  and 
creeping  waves  are  examples  of  some  of  the  possible  scattering  and  reflection 
mechanisms  [Ref.  l;p.  301].  The  "characteristic  impedance"  equals  the  density  of  a 
material  multiplied  by  the  speed  of  sound  in  the  material.  It  is  one  measure  used  to 
calculate  reflection  factors  for  pressure  waves.  Other  properties  that  affect  the 
scattering  properties  of  an  object  include  frequency,  shape  of  the  object,  and  pulse 
length. 

Ever  since  the  development  of  the  submarine  as  an  important  force  in  naval 
warfare,  much  effort  has  been  directed  to  improving  their  detection.  The  same 
holds  true  for  the  submarine  detecting  of  surface  ships  or  other  submarines.  The 
knowledge  required  for  detection  include  knowing  the  source  level  of  the  target, 
propagation  loss,  noise  level ,  and  lastly  target  strength,  the  subject  of  this  paper. 
These  parameters  are  used  in  the  well-known  sonar  equations  [Ref.  2:p.  411]. 
Target  strength,  TS,  as  used  in  the  sonar  equations  is  defined  as; 

TS=101ogi  (1) 

where  Ir  ia  the  reflected  intensity  one  meter  from  the  acoustic  center  of  the  target 
and  lo  is  the  incident  intensity  at  the  target  without  the  presence  of  the  target.  This 
equation  is  general,  giving  no  stipulation  as  to  the  direction  of  the  reflected 
intensity.  Typically,  signals  in  active  sonar  are  sent  and  recevied  by  the  same 
platform,  thus  the  most  common  geometry  and  ,therefore,  the  most  studied,  is  the 
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backscattered  monostatic  case  where  the  reflected  intensity  is  measured  in  the 
direction  back  toward  the  source.  Measurement  along  any  other  direction 
constitutes  a  bistatic  geometry.  Relatively  little  effort  has  been  given  to  the  study  of 
bistatic  target  strength  until  recently  when  the  bistatic  geometry  was  recognized  as  a 
potential  method  for  improving  underwater  detection  of  submerged  targets. 
Figures  1  and  2  illustrate  the  monostatic  and  bistatic  geometries. 


Source 

and 

Receiver 


Figure  1.  Backscattered  Monostatic  Geometry 


Source 


Figure  2.  Bistatic  Geometry 
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Some  work  has  been  done  in  the  radar  community  on  estimating  the  bistatic 
radar  cross- section  of  objects.  Crispin  and  Siegel  [Ref.  3:p.  156]  present  a  theorem 
for  the  approximation  of  radar  bistatic  cross-section: 

In  the  limit  of  vanishing  wavelength,  the  bistatic  cross-section  for 
transmitter  direction  ^  and  receiver  n  is  equal  to  the  monostatic  cross-section 
for  the  transmitter-receiver  direction  k  +  n  with  tc  /  ft  for  bodies  that  are 
sufficiently  smooth. 


(Figure  3  illustrates  the  geometry  as  stated  in  the  theorem.) 


Figure  3.  Radar  Bistatic  Cross-Section  Theorem  Geometry 
Their  bistatic  theorem  has  its  limitations  wHen  applied  to  complex  shapes;  The 
bistatic  cross-section  can  be  greater  or  less  than  the  monostatic  case  due  to 
structurally-related  interferences.  Another  limitation  is  the  size  of  the  target  with 
respect  to  the  incident  wavelength.  Spheres  were  analyzed  using  various  ka  values 
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where  where  a  is  the  radius  of  the  sphere  and  k  is  the  wave  number.  It  was  found 
that  the  theorem  held  more  closely  for  higher  ka  values,  but  in  all  cases  the  theorem 
sharply  deviated  from  experiment  as  the  bistatic  angle  approached  180  degrees.  As 
an  illustrative  example,  they  presented  data  on  a  jet  aircraft  and  a  ballistic-missile 
shape  [Ref.  3:p,  174-177].  Data  were  in  fairly  good  agreement  with  their  theorem 
as  long  as  the  bistatic  angle  did  not  get  too  large. 

In  transition  from  the  radar  community  to  underwater  acoustics,  Urick  [Ref. 
l:p.  302]  implied  that  the  bistatic  radar  cross-section  theorem  might  apply  in  a 
underwater  condition,  in  which  case  the  bistatic  target  strength  was  equal  to  the 
monostatic  target  strength  taken  at  the  bisector  of  the  bistatic  angle.  Urick 
emphasized  there  was  a  lack  of  data  to  prove  or  disprove  such  a  theory.  It  would  be 
very  convenient  if  such  a  theorem  existed  since  there  is  plentiful  target  strength  data 
for  monostatic  geometries.  Whether  or  not  the  high  frequency  limitations  of  such  a 
theorem  would  preclude  its  tactical  use  is  another  question. 

The  interest  in  bistatic  target  strength  has  been  generated  by  the  possibility  of 
increasing  detection  of  submerged  and  surface  units  using  this  technique.  For 
example,  one  can  easily  visualize  having  a  surface  ship  using  its  active  sonar  while  a 
submerged  submarine,  using  passive  sonar,  listens  for  reflections  off  enemy  forces. 

There  are  several  advantages  and  disadvantages  to  bistatic  employment.  Several 
ideas  now  discussed  were  taken  from  Harvey  [Ref.  4]  who  gave  a  thorough 
discussion  on  using  bistatic  employment  as  a  tactical  tool  in  naval  warfare.  TTie 
major  advantages  include  the  ability  to  maintain  a  track  of  the  enemy  force  using 
passive  sonar  and  the  increased  amount  of  data  available  to  detect  and  classify 
targets.  All  units  besides  the  active  unit  can  be  passive,  thereby  hampering  own- 
ship  detection.  Detecting  bistatic  reflections,  as  well  as  backscattered  reflections. 
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from  a  target  increases  the  data  available  for  detection  and  classification  of  the 
target.  The  greatest  difficulty  using  a  bistatic  geometry  is  the  high  positional 
accuracy  required  between  the  pinger  and  the  listener.  This  requires  prior  planning 
for  use  with  surface  ships  and  submarines  since  they  do  not  normally  have 
continuous  open  communication  with  each  other.  Another  difficulty  is  the 
distinguishing  of  backscattered  reflections  from  bistatic  reflections.  It  may  turn  out 
that  the  advantages  do  not  outweigh  the  disadvantages  and  the  whole  idea  will  prove 
infeasible. 
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A.  INTRODUCTION 


This  chapter  presents  some  definitions  and  background  knowledge  essential  in 
preparation  for  experimentation.  Some  concepts  are  rather  simple  but  are  included 
for  completeness. 

B.  DEFINITION 

The  definition  of  target  strength  has  already  been  given,  equation  (1). 
However,  there  are  several  aspects  of  the  definition  that  warrant  clarification. 
First,  all  measurements  are  taken  in  the  far  field  of  the  target.  The  far  field  is 
where  the  pressure  reflected  from  an  object  exhibits  a  1/r  relationship  [Ref.  2;p. 
1 78].  For  example,  pressure  measurements  taken  one  meter  from  the  hull  of  a  large 
ship  would  certainly  be  in  the  near  field,  so  intensity  is  measured  in  the  far  field  at 
ranges  of  hundreds  of  meters  and  then  extrapolated  back  to  one  meter  for  use  in  the 
target  strength  calculations. 

The  next  question  concerns  the  defintion  of  the  acoustic  center.  Very  rarely  are 
we  able  to  know  the  location  of  a  target's  acoustic  center.  The  acoustic  center  is  the 
point  from  which  sound  appears  to  radiate.  It  may  be  anywhere  around  the  target, 
and  may  lie  outside  the  shape  itself.  It  may  also  change  with  the  aspect  of  the  target 
relative  to  the  receiver.  Exact  location  of  the  target's  acoustic  center  is  ignored  in 
field  measurements,  and  the  geometrical  center  is  used.  When  measuring  pressure 
at  hundreds  of  meters,  a  few  meters  error  in  determination  accounted  for  by  the 
acoustic  center  will  not  produce  a  large  discrepancy. 


C.  SPHERICAL  WAVES 

For  a  distant  source,  to  a  good  approximation,  the  wave  will  appear  locally  as  a 
plane  wave  as  it  encounters  the  target.  Consider  an  omnidirectional  point  source: 
Spherically  diverging  waves  are  emitted  and  the  pressure  amplitude  and  phase  is 
dependent  only  on  the  distance  r  from  the  acoustic  center  of  the  source.  An 
outwardly  travelling  spherical  wave  can  be  represented  by 


P(r)=Ae' 

r 


(2) 


A  =  pressure  amplitude  at  one  meter 
r  =  distance  from  point  source 
w  =  angular  frequency  of  signal 
k  =  wave  propagation  number 

When  the  radius  of  curvature  of  the  wave  front  is  large  compared  to  the  biggest 
dimension  of  the  target,  the  spherical  wave  can  be  treated  as  a  plane  wave. 

Plane  waves  are  assuired  when  the  din^nsions  of  a  target  are  small  compared  to 
the  source-target  distance.  One  of  the  goals  of  small-scale  experiments  is  the 
application  of  the  results  to  oceanic  situations,  therefore  there  is  a  concern  about  the 
approximation  of  plane  waves  ensonifying  a  target  in  a  small  laboratory  tank. 

The  geometrical  constraints  of  the  water  tank  limits  target-source-receiver 
distances.  With  small  distances  between  source  and  receiver  and  big  targets,  a  plane 
wave  may  not  be  a  good  approximation.  Another  factor  that  influences  the  plane 
wave  approximation  is  the  directionality  of  the  source.  It  is  difficult  to  achieve 
omnidirectional  behavior  from  a  finite  source  and  it  is  impossible  if  the  frequency 
is  high.  Beam  patterns  are  formed  with  energy  emitted  in  preferred  directions.  If 
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the  target  subtends  an  angle  greater  than  the  beamwidth,  the  pressure  amplitude  will 
not  be  the  same  for  all  points  on  the  target. 

Plane  waves  also  have  properties  that  make  target  strength  calculations  simple 
to  handle.  The  intensity  can  be  equated  to  the  root-mean-square  pressure  amplitude 
squared  divided  by  the characterisdc  impedance  [Ref.  2;p.  111].  Since  hydrophones 
measure  pressure  via  a  sensitivity  factor,  the  output  voltage  they  produce  can  be 
directly  related  to  intensity. 

D.  WAVEFORM  EFFECTS 

There  are  several  factors  to  consider  in  the  measurement  and  collection  of 
reflected  waveforms  from  objects.  Some  are  listed  below  and  each  will  cause  the 
reflected  waveform  to  differ  from  the  incident  waveform. 

-  target  composition 

-  target  shape 

-  pulse  length 

-  resonance  effects 

For  the  reflected  waveform  to  be  identical  of  the  incident  waveform,  a  planar 
rigid  boundary  must  be  encountered.  For  more  complicated  targets,  some  of  the 
incident  energy  will  be  transmitted  into  the  body  and  eventually  reradiated  and 
different  areas  of  the  target  may  reflect  in  such  a  way  to  cause  mutual  interference. 
Many  objects  have  resonances,  at  which  the  target  will  store  some  of  the  incident 
energy  and  reradiate  it  at  a  later  time  causing  interference  effects.  Theoretical  and 
experimental  results  show  that  for  short  pulse  lengths,  reflections  are  nearly 
independent  of  target  composition  and  are  almost  identical  to  the  incident  pulse 
[Ref.  5,6,7]  while  for  longer  pulses  the  reflected  waveform  changes  drastically 
depending  on  the  composition  and  shape  of  the  target. 
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E.  EXPERIMENTAL  MEASUREMENT 


While  conducting  the  literature  search,  it  was  interesting  to  note  how  the 
pressure  amplitude  of  the  reflected  wave  was  measured.  As  stated  before, 
hydrophones  give  a  voltage  that  is  representative  of  pressure  and  this  can  be  equated 
to  intensity.  The  problem  lies  in  the  fact  that  the  reflected  waveform,  in  most  cases, 
is  not  identical  to  that  of  the  incident  signal.  The  experimenter  has  to  specify  what  is 
measured  and  what  kind  of  pulse  is  being  used.  Response  can  be  catagorized  as 
transient  or  steady-state.  Steady-state  implies  that  the  whole  target  has  been 
ensonified  by  the  incident  pulse  allowing  for  all  portions  of  the  target  to  contribute 
to  the  reflection  and  for  any  resonances  to  become  fully  excited.  In  this  case,  the 
return  will  be  a  combination  of  all  possible  reflection  paths  and  modes.  A  transient 
response  simply  is  not  in  the  steady-state.  To  give  an  idea  of  measurement  used,  the 
following  are  quoted  from  two  published  articles: 

...five  cycles  of  the  sine  wave  located  at  a  position  in  the  pulse  judged  to 
be  a  true  representation  of  the  steady  state.  The  average  of  the  ten  peak  values 
so  obtained  gave  the  peak  pressure  amplitude  [Ref.  7] 


After  amplifrcation,  the  amplitude  of  the  back- scattered  signal  is 
measured  with  a  peak  detector  [Ref.  8]. 

The  main  point  is  that  for  reflected  waveforms,  not  similiar  in  shape  to  the 
incident  waveform,  some  subjective  determination  is  required.  Estimation  of  the 
steady-state  must  be  made  or  some  definite  parameter  has  to  be  selected  such  as  peak 
amplitude  or  total  energy.  In  any  case,  an  unique  way  to  convert  the  amplitude 
measurements  to  an  intensity  level  is  not  achievable.  Add  to  this  the  fact  that  the 
reflected  waveform  will  be  different  for  different  target  orientations  and 
frequencies,  and  the  problem  of  translating  measurements  into  an  intensity  level  to 
calculate  a  target  strength  becomes  complicated. 


There  presently  exists  a  Target  Strength  Standardization  Program  [Ref.  9] 
which  was  assigned  the  task  of  finding  some  standard  measure  of  target  strength 
independent  of  signal  type.  The  general  consensus  from  the  program  was  to  use  the 
"Average  Integrated  Target  Strength"  which  is  the  average  time-integrated  energy. 
Furthermore,  an  Echo  Metrics  Subcommittee  [Ref.  9]  was  formed  to  find  signal- 
processor  dependent  metrics  for  the  measurement  of  target  strength.  Theoretically, 
the  incident  energy  encountered  by  the  object  would  be  reflected  or  absorbed  and 
eventually  reradiated  at  some  later  time  providing  no  energy  was  stored  or  lost  in 
the  object.  Therefore,  it  would  seem  that  the  average  of  the  integrated  energy 
would  be  an  independent  metric  for  all  types  of  signals  providing  the  integration 
time  is  extended  to  include  all  the  return  signal.  Some  standardization  of  target 
strength  is  being  investigated  so  that  data  from  independent  experiments  can  be 
compared  and  integrated. 
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m.  DATA  ACQUISITION  SYSTEM 

A.  BACKGROUND  AND  CONCEPT 

The  experiments  were  conducted  in  a  water  tank  located  in  Room  019  in 
Spanagel  Hall  at  the  Naval  Postgraduate  School.  The  tank  was  aj^oximately  1.S  m 
wide,  8  m  long,  and  2  m  deep  and  lined  with  anechoic  rubber  material.  For  the 
targets  used,  wall  reflections  were  of  the  same  magnitude  as  the  target  reflections. 
Usually,  it  was  easier  to  use  pulsed  signals  and  time  discriminate  for  the  target 
reflection.  However,  the  small  dimensions  of  the  tank  did  not  provide  the 
environment  to  do  this  effectively  for  all  geometries.  To  overcome  this  difficulty,  a 
subtraction  technique  was  used  employing  a  digital  storage  device.  The  basic 
concept  was  to  subtract  the  signal  received  with  the  target  immersed  in  the  tank 
from  the  signal  received  with  the  target  absent.  The  procedure  was  very  simple  and 
allowed  for  any  target-source-receiver  geometry. 

B.  DEVELOPMENT 

A  considerable  amount  of  time  went  into  the  selection  of  the  digital  sampling 
and  storage  device,  the  source,  and  the  receiver.  The  Hrst  question  that  had  to  be 
addressed  was  the  range  of  frequencies  to  be  used.  The  frequency  range  of  100-200 
kHz  was  initially  selected  because  it  gave  wavelength-to-target  dimension  ratios  that 
match  oceanic  scenarios.  The  critical  specifications  for  the  digital  storage  device 
are  the  sampling  rate,  dynamic  range  and  jitter.  Dynamic  range  was  not  a  problem 
since  most  digital  oscilloscopes  have  the  same  voltage  range  as  a  typical  non-digital 
oscilloscope.  The  sampling  rate  must  be  fast  enough  to  provide  a  good 
representation  of  the  signal,  and  assuming  a  sinusoidal  signal  about  ten  samples  per 
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cycle  would  give  a  good  trace.  For  example,  for  a  100  kHz  signal  the  sampling  rate 
would  have  to  be  1  MHz.  Another  major  concern  was  the  amount  of  jitter  in  the 
sampling  device.  Jitter  describes  the  fluctuation  in  the  triggering  times  relative  to 
the  sampled  signal.  For  a  system  using  a  subtraction  process,  it  is  perhaps  the  most 
critical  specification  because  there  has  to  exist  a  common  time  reference  between 
collected  waveforms  to  get  accurate  results. 

The  selection  of  the  digital  oscilloscope  was  based  on  the  availability  of 
equipment.  Two  digital  oscilloscopes  were  evaluated.  The  Nicolet  3091  was 
analyzed  first.  Its  fastest  sampling  rate  was  1  MHz  and  had  a  storage  capability  of 
4000  points.  The  Nicolet  3091  had  a  RS-232  serial  interface  and  was  easy  to  use, 
however  it  was  not  acceptable  due  to  excessive  jitter  at  the  frequencies  of  interest. 
The  jitter  was  easily  seen  by  applying  a  pulsed  signal  and  observing  the  traces  on  the 
scope.  The  digital  oscilloscope  finally  used  was  the  Tektronix  2430  (Tek  2430).  It 
is  a  very  capable  oscilloscope  and  has  many  nice  features  making  it  suitable  for  data 
acquisition.  The  most  important  feature  that  it  had  was  an  averaging  mode 
selectable  from  2  to  256  waveforms.  This  helped  eliminate  random  noise  and 
minimized  any  fluctuations  the  electronic  equipment  possessed.  The  amount  of 
jitter  was  dependent  of  the  time  scale  selected.  The  fastest  sampling  rate  of  the  Tek 
2430  was  100  MHz  which  could  be  used  to  sample  a  10  MHz  signal  using  ten  points 
per  cycle.  The  jitter  specification  given  in  the  operating  manual,  for  the  time- 
delayed  mode,  was  .08  x  (time  scale).  The  jitter  was  observed  and  was  signiflcant 
for  the  time  scale  desired,  however  the  saving  factor  was  the  time  averaging  mode 
which  ultimately  averaged  out  the  jitter  and  provided  a  sufficiently  constant  time 
reference.  It  was  very  easy  to  observe  this  on  the  oscilloscope  because  of  a 
triggering  marker  on  the  signal.  Jitter  becomes  less  of  a  problem  for  any  digital 
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oscilloscope  when  the  number  of  samples  per  cycle  of  signal  increases.  However, 
the  problem  then  arises  from  the  limited  amount  of  storage  in  the  oscilloscope.  A 
smaller  number  of  cycles  can  be  stored  when  using  a  faster  sampling  rate.  The  Tex 
2430  has  a  storage  capabiltiy  of  1024  points,  much  less  that  the  Nicolet  3091; 
however  it  also  has  a  time-delay  feature  which  allowed  for  selective  data  storage 
(The  Nicolet  3091  did  not  have  any  time-delay  feature).  Other  features  of  the  Tek 
2430  include  cursor  control  for  voltage  and  time  measurements,  simple  output  to  a 
Hewlett  Packard  Think-Jet  printer  or  pen  plotter,  and  a  GPIB  interface.  It  is  a  very 
versatile  piece  of  equipment  and  was  finally  selected  for  use  in  the  data  acquisition 
system. 

The  selection  of  the  source  and  receiving  hydrophone  was  also  limited  to  in- 
house  equipment.  For  the  frequencies  desired,  the  USRD  type  F41  transducer  was 
the  appropriate  choice  with  a  nominal  frequency  range  of  15-150  kHz  [Ref.  lOq).  5]. 
It  is  a  piston  type  piezeoelectric  transducer  with  a  face  diameter  of  10  cm.  The 
choice  of  the  receiver  was  a  Celesco  LC-IO  hydrophone.  It  has  a  wide  nominal 
frequency  band,  is  small  and  light  making  it  easy  to  move  around,  and  has  good 
sensitivity. 

C.  SYSTEM  SET-UP 

1.  Electronic  Equipment 

The  F41  source  had  to  be  driven  near  200  V  rms  for  target  reflections  to  be 
seen.  An  Instruments,  Inc.  Model  L2  power  amplifier  was  used  to  drive  the  source 
and  its  selection  was  again  based  on  what  was  available  in  the  laboratory.  It  allowed 
for  impedance  matching  and  has  a  voltage-divided  10  mV/V  port  to  monitor  output 
voltage.  The  output  of  the  receiver  was  passed  through  an  Ithaco  Model  1201  Low 
Noise  preamplifier  which  has  bandwidth  setting  capability.  This  was  very 


beneficiai  because  there  was  a  lot  of  60  Hz  noise  which  had  to  be  eliminated  by 
bandwidth  Altering.  The  signal  out  of  the  preamplifier  was  sent  to  a  Kikusui 
COS5060A  oscilloscope  and  the  Tektronix  2430  digital  oscilloscope.  Having  a  non¬ 
digital  oscilloscope  made  observing  and  measurement  of  the  waveforms  easier  to 
perform  on  the  Tek  2430.  A  Hewlett  Packard  3314A  function  generator  was  used 
to  provide  the  source  signal.  It  also  provided  the  triggering  signal  which  was  sent  to 
the  external  channels  of  the  oscilloscopes.  In  the  pulse  mode,  the  HP  3314A  sends 
out  a  3  V  transistor-transistor  logic  signal  via  its  trigger  port.  The  pulse  repetition 
rate  of  10  ms  was  sufficient  to  allow  reflections  to  diminish  before  the  next  pulse.  A 
Hewlett  Packard  9000  Series  300  desk  top  compute  was  used  as  the  data  collection 
and  processing  unit.  It  was  connected  to  the  Tek  2430  via  the  HPIB  bus.  The 
transfer  of  waveforms  stored  on  the  Tek  2430  was  controlled  by  software 
programs.  Data  was  ASCII  encoded  from  the  Tek  2430,  transferred,  decoded  into 
actual  voltages,  and  then  stored  onto  disks  in  binary  data  form.  A  schematic 
diagram  of  the  electronic  system  is  shown  in  Figure  4. 

2.  PhysitaL  Sstrttp 

To  measure  all  possible  bistatic  angles,  the  scattering  pattern  would  have  to 
obtained  in  three  dimensions.  This  was  not  done  in  this  work.  Bistatic  data  were 
limited  to  the  horizontal  plane  determined  by  the  source,  receiver,  and  target,  with 
the  axis  of  the  target  in  the  same  plane. 

The  source,  the  F41  transducer,  was  mounted  on  a  3/4-inch  aluminum  rod 
attached  to  a  rectangular  sufqxirt  plate.  The  rod  was  fixed  to  the  plate  by  means  of  a 
friction  Atting  secured  by  a  hose  clamp.  Target,  source,  and  receiver  were 
submerged  one  meter  below  the  surface  of  the  water  and  the  rods  which  supported 
them  were  marked  for  this  depth.  The  source  with  its  support  rested  on  another 
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HP  3314A 
Function  Generator 


support  built  from  aluminum  L-beams  that  spanned  the  width  of  the  tank.  Lead 
bricks  were  used  to  hold  the  source  once  it  was  set  in  place  and  a  level  was  used  to 
ensure  the  source  face  was  vertical.  The  receiver  was  mounted  onto  a  sliding  track 
specially  built  for  the  tanks.  Adjusting  screws  were  added  for  ease  of  positioning 
the  receiver  vertically.  The  inidal  intent  was  to  design  a  rotating  boom  to  hold  the 
receiver  and  support  the  target  from  its  center.  A  device  was  constructed,  but  too 
much  receiver  movement  occurred  between  the  loading  and  unloading  of  targets. 
Though  the  boom  was  never  used,  the  device  was  still  used  for  the  suspension  of 
targets  and  had  a  mechanical  angle  indicator  on  its  top  for  angle  measurement. 
Spherical  targets  were  suspended  by  a  thin  cable  and  the  cylindrical  target  by  a  3/8- 
inch  rod. 

The  direction  of  the  source  was  set  by  aligning  the  source,  receiver,  and 
target  along  a  straight  line  with  the  receiver  between  the  source  and  target.  Source, 
receiver,  and  target  depth  were  one  meter.  Vertical  alignment  was  ensured  by 
keeping  the  source  and  receiver  vertical  by  use  of  a  level.  The  frequency  was  set  to 
200  kHz  to  establish  a  narrow  major  beam  and  the  source  was  rotated  until 
maximum  signal  was  seen  by  the  receiver.  This  ensured  the  target  was  being  totally 
ensonified  by  the  major  lobe  of  the  source. 

Once  the  source  direction  was  set,  the  receiver  was  moved  to  the  desired 
position.  The  angle  between  the  source,  target,  and  receiver  was  determined  by 
setting  the  receiver's  track  perpendicular  to  the  source-target  axis  and  then  using 
trigonometry  relationships.  After  the  desired  geometry  was  set,  care  was  taken  to 
ensure  the  target,  source,  and  receiver  were  stationary.  It  was  critical  that  all 
elements  be  stationary  when  waveforms  were  collected. 
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3. 


The  method  of  data  collection  and  storage  followed  the  sequence; 

-  set  up  desired  source-receiver-target  geometry 

-  warmup  all  electronic  equipment  >1  hr 

-  ensure  connections  as  illustrated  in  Figure  5 

-  load  "Tex"  onto  computer 

-  set  up  the  Textronix  2430  oscilloscope 

-  B  trigger  mode  to  "runs  after" 

-  A  trigger 

-  mode:  nonmal 

-  coupling:  AC 

-  source:  external  1 

-  slope:  positive 

-  level:  800  mV 

-  trigger  position:  1/8 

-  data  encoding:  ASCII 

-  GPIB  address:  4 

'  vertical  mode:  channel  2 
'  select  time  delay  to  obtain  desired  time  window 

-  set  acquisition  mode  to  256  averaging 

-  adjust  channel  2  voltage  scale/Ithaco  gain  appropriately 

for  the  signal 

-  for  each  new  signal,  press  the  acquire  button 

-  wait  sufficient  time  for  averaging  to  take  place 

-  press  the  save  button 

-  run  the  "Tex"  program 

The  "Tex"  program  controlled  the  transfer  and  storage  of  data  from  the 
digital  oscilloscope.  Inputs  from  the  user  were  file  names  and  ^le  descriptions. 
The  transfer  process,  including  writing  onto  the  disk,  took  af^roximately  10  s  per 
waveform.  A  listing  of  the  program  is  in  Aj^ndix  A. 
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The  computer  program  "Plot"  was  used  to  evaluate  the  stored  data.  It 
took  two  files  and  subtracted  them,  and  then  plotted  the  results  on  the  CRT 
terminal.  Hard  copy  output  could  be  obtained  by  porting  to  a  pen  plotter  or  Think- 
Jet  printer.  Input  from  the  user  required  file  names  and  where  printout  was  to 
occur.  The  plotting  routine  was  written  so  that  the  plot's  parameters  were  easily 
changeable.  A  listing  is  enclosed  in  Appendix  B. 

The  method  of  data  collection  was  a  slow  and  somewhat  tedious  process 
since  all  movement  of  target  and  receiver  ware  manually  controlled.  Data 
readings  were  also  recorded  manually. 

4.  Validation 

The  question  of  whether  the  data  acquisition  system  could  correctly 
extract  a  target  reflection  immersed  within  other  signals  had  to  be  answered.  The 
very  simple  procedure  of  subtracting  background  noise  works  only  if  there  is  a 
common  time  reference  for  "target  in"  and  "target  out"  waveforms.  For  ease  of 
discussion,  "background  noise"  is  used  for  any  signal,  random  or  coherent,  that  is 
not  the  desired  signal  and  corresponds  to  the  "target  out"  condition. 

The  procedure  to  validate  the  system  began  with  obtaining  a  geometry 
which  allowed  a  direct  pulse  and  reflected  pulse  to  be  completely  separated  by 
time  discrimation  with  no  interfering  wall  reflections.  Figure  5  is  such  an  example. 
The  larger  pulse  was  the  direct  path  pulse  and  the  second  pulse  is  the  pulse 
reflected  from  a  styrofoam  sphere.  The  frequency  was  66.5  kHz  and  the  pulse 
length  was  20  cycles.  The  reflected  pulse  is  shown  enlarged  in  Figure  6.  The 
peak-to-peak  voltage  of  the  reflected  pulse  was  measured  to  be  68  ±  2  mV  using 
the  cursor  controls  of  the  Tek  2430. 
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The  pulse  length  was  then  increased  to  120  cycles  so  the  direct  pulse  would 
be  interfering  with  the  reflected  pulse.  No  other  changes  were  made.  The 
interference  of  the  two  signals  is  shown  in  Figure  7  and  was  stored  in  file  A.  Next, 
the  sphere  was  removed  and  the  background  noise  sampled  and  recorded  in  file  B. 
The  background  noise  in  this  case  was  mostly  the  direct  pulse.  The  last  step  was  to 
subtract  file  B  from  file  A  to  get  the  reflected  signal  from  the  sphere.  The  result  is 
shown  in  figure  8.  The  waveform  begins  at  the  same  time  in  comparison  with 
Figure  6  and  the  only  difference  is  that  the  pulse  length  was  longer  as  expected  and 
the  curve  slightly  more  choppy  in  form.  Measured  voltage  was  73  ±  4  mV.  The 
larger  uncertainty  arose  due  to  the  combination  of  small  jitter  and  amplitude 
fluctuations  of  the  HP  33 14A.  The  uncertainty  was  estimated  beforehand  by  taking 
identical  signals  directly  from  the  function  generator  to  the  Tek  2430,  recording 
them  in  separate  files,  and  subtracting  them  and  observing  the  results.  It  was  found 
that  peak-to-peak  voltage  uncertainty  was  approximately  four  digitizing  levels 
regardless  of  the  voltage  scale  used  (a  digitizing  level  equals  the  voltage  per  division 
divided  by  50  for  the  Tek  2430).  Given  that  this  error  would  always  be  present, 
there  was  good  agreement  between  Figure  6  and  Figure  8,  and  the  method  of 
obtaining  data  was  considered  validated. 

As  with  any  signal  processing  system,  a  high  signal-to-noise  ratio  is 
desired.  For  any  digital  subtraction  routine,  the  desired  signal  has  to  be  large 
enough  relative  to  the  noise  so  it  will  not  disappear  in  the  digitizing  level  of  the 
sampling  device.  For  example,  if  noise  plus  the  signal  were  10  volts  using  512 
digitizing  levels  or  bins,  each  bin  represents  19,5  mV.  If  the  signal  were  small,  i.e., 
1  mV,  the  signal  level  would  be  much  smaller  that  the  digitizing  level  and  would  be 
lost  in  the  subtraction  routine. 


19 


Keeping  the  source  and  receiver  positions  fixed  relative  to  each  other  was 
critical.  Small  movement  between  the  two  greatly  affected  the  phase  of  the  received 
signal  and  had  a  large  effect  on  the  results.  This  effect  was  verified  in  preliminary 
work.  The  initial  intent  was  to  place  the  receiver  and  support  the  target  on  the  same 
platform.  However,  the  removal  of  the  target  for  the  background  noise 
measurement  caused  flexing  in  the  platform  and  produced  source-to-receiver 
movement  giving  incorrect  results  in  the  subtraction  routine.  Figure  9  illustrates 
the  effect  of  source-receiver  movement  on  "target-in"  and  "target-out"  waveforms. 
The  solid  trace  represents  the  target-in  waveform  and  the  dashed  line  represented 
the  target-out  waveform.  The  reflection  began  at  3.5  usee.  The  dashed  and  dotted 
line  is  the  the  result  from  the  subtraction  of  the  two  waveforms  and  it  was  clear  that 
in  this  case,  it  did  not  give  the  correct  results.  The  region  between  0  and  35  usee 
should  have  cancelled  to  zero.  The  non-zero  result  indicates  a  phase  change  caused 
by  relative  source-receiver  movement.  The  phase  change  subsequently  affected  the 
amplitude  of  the  reflected  signal  from  35  to  100  usee.  It  was  easy  to  see  that  any 
phase  change  brought  about  by  small  source-receiver  movement  would  cause 
erroneous  results.  To  eliminate  this  problem,  a  separate  platform  was  subsequently 
built  and  used  for  the  receiver. 

The  output  of  the  power  amplifier  was  also  checked  in  preliminary  work 
for  fluctuations.  An  arbitrary  source-receiver  geometry  was  selected  without  any 
target  in  the  tank.  Several  waveforms  of  the  direct  pulse  were  stored  at  various 
times  after  initial  start  up.  Amplitude  fluctuations  existed  shorty  after  start-up. 
Ideally,  successive  waveforms  subtracted  from  each  other  would  cancel  out  to  give  a 
zero  response.  If  amplitude  fluctuations  existed  a  false  signal  would  occur  and  this 
in  fact  did  occur  after  five  minutes  of  warm-up  as  Figure  10  illustrates.  Another 
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test  taken  at  15  minutes  after  start  up  produced  better  but  still  unsatisfactory  results. 
After  one  hour  of  warm-up  the  amplitude  fluctuations  subsided  and  a  zero  respose 
was  obtained. 


Pe  I  at  t  ve  T i me 


Figure  10.  Waveform  Amplitude  Fluctuation  Effect  After  10  Minutes 

of  Warm-up 


IV. 


TTie  measurement  of  the  far  field  of  several  objects  was  conducted.  An  equation 
that  gives  the  approximate  beginning  of  the  far  field  is  given  by  [Ref.  2:p.  188]: 


L  =  the  longest  dimension  of  object 
A.  =  wavelength 

It  is  always  best  to  determine  the  far  field  by  measuring  pressure  amplitude  as  a 
function  of  distance  from  the  object.  The  equation  for  a  diverging  spherical  wave, 
as  described  by  equation  2,  can  be  manipulated  by  taking  the  log  of  both  sides  to 
transform  it  into  a  linear  equation  with  a  slope  of  minus  one.  The  far  field  can  be 
determined  experimentally  when  the  observed  slope  of  the  curve  starts  to  deviate 
from  minus  one. 

The  far  field  of  the  F41  was  determined  to  begin  at  approximately  100  cm  and 
70  cm  for  176  kHz  and  70  kHz,  respectively.  The  far  field  for  the  sound  reflected 
from  the  styrofoam  cylinder  positioned  broadside  to  the  source  was  determined  to 
begin  at  40  cm  for  70  kHz.  The  three-inch  diameter  styrofoam  sphere's  far  field 
began  at  less  than  20  cm  at  70  kHz.  Figures  1 1  through  14  are  graphs  of  the  data 
used  to  determine  their  far  fields.  Lines  of  slope  equal  to  minus  one  are  drawn  on 
the  graphs. 
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V.  SPHERES 


A.  INTRODUCTION 

To  achieve  an  understanding  of  the  complexities  involved  in  analyzing 
scattering  effects  of  various  objects,  it  was  instructive  to  begin  with  a  very  simple 
and  symmetrical  object,  a  sphere.  Several  solid  metal  spheres  were  used.  Spheres 
were  aluminum  and  brass,  but,  the  exact  composition  was  not  known  or  investigated 
since  they  were  to  serve  only  as  preliminary  tools  in  studying  the  general  properties 
of  waveform  scattering. 

Since  metal  spheres  are  not  perfectly  rigid,  they  are  capable  of  absorbing 
energy  and  subsequently  reradiating  it.  Also,  they  possessed  resonances.  Many 
experiments  in  the  past  30  years  have  been  conducted  with  spheres  and  long 
cylinders,  and  the  Resonance  Scattering  Theorem  (RST)  has  evolved  from  this  work 
[Ref.  5,6,7,8,11,12,13,14].  This  theorem  predicts  scattering  patterns  which  entail 
both  the  rigid-like  properties  and  the  resonance  properties  of  an  object.  It  was  not 
the  intent  to  examine  this  theorem  in  detail  but  only  to  observe  some  of  its 
predictions  experimentally. 

B.  MONOSTATIC  FORM  FUNCTION 

The  goal  of  replicating  some  of  the  data  in  the  literature  was  sought  to  establish 
confidence  in  the  experimentation  process.  Neubauer,  Vogt,  and  Dragonette  [Ref. 
7]  conducted  an  analysis  on  elastic  spheres  using  steady-state  signals  and  observed 
the  backscattered  signal  form  function  from  an  aluminum  sphere  as  a  function  of  ka. 
The  reflection  form  function  is  defined  as: 
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(4) 


f  =  ^  Ps 
a  Po 

Ps  =  scattered  pressure  amplitude  measured  at  r 
Po  =  incident  pressure  amplitude  at  sphere  center 
a  =  sphere  radius 

Their  experimental  data  were  plotted  and  compared  against  theoretical  results 
and  good  agreement  occurred. 

Part  of  Neubauer's  experiment  was  repeated.  A  four-inch  diameter  aluminum 
sphere  was  ensonified  by  50  cycle  pulses  from  100  to  130  kHz,  and  the  backscattered 
pressure  amplitude  was  measured.  Time  discrimination  of  the  reflection  was  used. 
The  amplitude  of  the  reflected  pulse  was  taken  at  what  was  judged  to  be  the  steady- 
state  part  of  the  signal.  The  receiver  was  1.60  m  from  the  center  of  the  sphere, 
placing  the  receiver  well  within  the  far  field  of  the  sphere. 

The  wavelength  used  in  the  calcuation  of  ka  was  dependent  on  the  speed  of 
sound  in  the  water,  which  in  turn  was  affected  by  temperature,  pressure,  and 
salinity.  For  the  fresh  water  tanks,  salinity  was  ignored.  The  speed  of  sound  was 
based  on  an  empirical  formula  that  applied  to  distilled  water  [Ref.  l:p.  107]. 
Although  the  water  was  normal  tap  water,  the  distilled  water  formula  was  a  close 
enough  approximation  [Ref.  15].  Water  temperature  of  the  tank  was  a  consistent  17 
degrees  Celsius  throughout  the  entire  experimentation  phase.  At  this  temperature 
and  normal  atmospheric  pressure,  the  speed  of  sound  was  calculated  to  be  1473  m/s. 
Experimentally,  two  LC-lO's  were  placed  on  a  track  180  cm  apart  and  the 
successive  arrival  times  of  the  pulse  were  measured  via  the  Tek  2430  with  the 
cursor  controls.  The  speed  of  sound  thus  was  measured  to  be  1480  ±  20  m/s,  placing 
1473  m/s  well  within  the  experimental  uncertainty,  and  this  value  was  used  as  the 
speed  of  sound  in  water. 
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The  shadowing  effect  of  the  LC-10  in  front  of  the  F41  source  was  taken  into 
account  by  placing  a  dummy  LC-10  at  that  position  and  recording  incident  pressure 
at  the  sphere  center  location.  Additionally,  since  the  LC-IO  was  not  perfectly 
omnidirectional,  the  support  rod  was  marked  to  maintain  a  constant  receiver  aspect 
with  the  target.  Placing  the  LC-10  in  front  of  the  F41  also  caused  additional 
reflections  to  occur.  The  direct  pulse  was  reflected  by  the  LC-10  back  to  the 
source's  face  and  again  reflected  there  back  to  the  LC-10.  More  "noise"  was  present 
when  measuring  the  monostatic  case  but  appropriate  location  of  the  receiver 
eliminated  any  interference  from  those  reflections. 

The  data  matched  Neubauer's  results  well,  and  the  form  function  as  a  function 
of  ka  is  plotted  in  Figure  15.  Data  from  Neubauer’s  experiment  is  also  plotted  in 
Figure  15.  The  only  difference  between  results  was  a  small  shift  in  ka  of 
approximately  0.2-0.3,  not  surprising  since  the  compostion  of  the  sphere  was 
unknown.  The  minima  corresponds  to  resonant  frequencies  of  the  sphere,  where 
the  specularly  reflected  wave  is  interfering  with  the  resonantly  radiated  wave  from 
within  the  body. 

C.  WAVEFORM  EFFECTS 

The  effect  of  a  sphere's  resonance  was  examined  in  more  detail.  A  four- inch 
diamter  aluminum  sphere  was  ensonified  by  a  30  cycle  pulse  at,  above,  and  below  a 
resonant  frequency  of  82.7  kHz.  At  resonance,  the  ringing  of  the  sphere  is  clearly 
seen  in  Figure  16  when  the  direct  pulse  ends.  Pulse  length  was  approximately  367 
|i,s  long.  The  sphere's  stored  energy  is  simply  being  reradiated.  Figures  17  and  18 
show  the  waveform  response  2  kHz  above  and  below  resonance.  They  are 
decisively  different  and  clearly  illustrate  the  resonance  effect. 
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The  effect  of  sphere  composition  on  reflected  waveforms  was  observed.  A 
four-inch  diameter  brass  sphere  was  ensonified  by  the  same  frequencies  as  the 
aluminum  sphere.  Brass  has  a  characteristic  impedance  approximately  2.5  times 
that  of  aluminum.  Figures  19,  20,  and  21  show  the  reflected  waveforms  having 
frequency  dependence  that  are  markedly  different  from  the  aluminum  sphere. 

The  effect  of  pulse  length  was  observed  using  a  4-inch  aluminum  sphere  at  and 
off  a  resonant  frequency.  Figure  22,  at  resonance,  and  Figure  23,  off  resonance, 
show  reflected  wavforms  using  a  short  pulse  of  2  cycles.  They  appear  very  similar 
in  shape.  Figures  24  and  25  show  the  results  for  a  pulse  length  of  5  cycles.  Again, 
the  results  of  the  two  are  similar.  When  pulse  length  reaches  15  cycles,  the  effects 
of  the  resonance  begins  to  be  clearly  observed,  as  Figure  26  and  27  illustrate. 

For  comparison,  the  same  procedure  above  was  performed  on  a  four-inch 
diameter  brass  sphere.  Results  are  illustrated  in  Figures  28  through  33.  The  brass 
sphere  did  not  show  any  resonance  effect  as  determined  in  the  preceding 
composition  experiment  for  these  frequencies,  so  its  results  were  very  similar  to  the 
aluminum  sphere’s  results  off  resonance.  Waveforms  were  practically  identical  for 
the  shorter  pulse  lengths  regardless  of  frequency  and  composition  implying  the 
reflections  are  purely  specular.  The  results  clearly  demonstrated  the  dependence  of 
pulse  length  on  reflected  waveforms. 

D.  SELECTION  OF  STYROFOAM 

Since  defining  target  strength  becomes  difficult  when  using  elastic  spheres, 
which  was  not  the  intent  of  this  work,  the  need  existed  to  do  one  of  two  things;  (1) 
work  at  a  frequency  that  gives  good  replication  of  the  input  pulse,  or  (2)  find  a 
target  that  gave  a  very  high  impedance  mismatch. 
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Figure  16.  Aluminum  Sphere  Reflection  at  Resonance,  82.7  kHz 
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Figure  17.  Aluminum  Sphere  Reflection,  84.7  kHz 
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Figure  18.  Aluminum  Sphere  Reflection,  80.7  kHz 
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Figure  20.  Brass  Sphere  Reflection,  82.7  kHz 
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Figure  21.  Brass  Sphere  Reflection,  84.7  kHz 


Figure  22.  Aluminum  Sphere  Reflection,  82.7  kHz 
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Figure  23.  Aluminum  Sphere  Reflection,  84.7  kHz 


Figure  25.  Aluminum  Sphere  Reflection,  84.7  kHz 
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Figure  26.  Aluminum  Sphere  Reflection,  82.7  kHz 
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Figure  27.  Aluminum  Sphere  Reflection,  84.7  kHz 
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Figure  30.  Brass  Sphere  ReHection,  82.7  kHz 
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Figure  31.  Brass  Sphere  Reflection,  84.7  kHz 
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Figure  32.  Brass  Sphere  Reflection,  82.7  kHz 
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Figure  33.  Brass  Sphere  Reflection,  84.7  kHz 
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The  most  important  factor  in  achieving  a  good  reflection  is  the  impedance 
mismatch  between  the  fluid  and  object.  Experimental  work  done  with  metal  spheres 
by  Hamptom  [Ref.  5]  shows  the  effect  of  compostion  using  several  metals.  A 
pressure  wave  in  water  will  see  little  resistance  to  a  hard  rubber  sphere  because  its 
characteristic  impedance  is  near  that  of  the  water.  An  object  that  gives  a  high 
impedance  mismatch  is  an  air  bubble.  The  impedance  ratio  of  water  to  air  is 
roughly  3(XX).  For  the  case  of  an  air  bubble,  the  boundary  condition  between  the 
water  and  air  is  termed  "pressure  release",  since  the  pressure  transmitted  into  the  air 
is  nearly  zero  [Ref.  2:p.  127].  The  difference  between  a  pressure-release  sphere  and 
a  rigid  sphere  is  best  seen  graphically  in  Figure  34  [adapted  from  Ref.  15:  p.  359]. 
In  the  ka  region  used  in  the  experiments,  both  the  rigid  and  pressure-release  form 
functions  will  be  nearly  identical  and  equal  to  one. 


Figure  34.  Backscattered  Form  Function  of  a  Rigid  and  Pressure 

Release  Sphere 

A  styrofoam  ball  was  used  to  approximate  an  air  bubble  not  near  resonance. 
The  reflected  waveforms  were  generally  identical  to  the  input  pulse  at  all 
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frequencies.  This  made  styrofoam  an  ideal  target  for  this  work.  Figures  35 
through  40  show  reflections  from  a  three-inch  diameter  styrofoam  sphere  using 
various  pulse  lengths  at  100  and  106  kHz. 

E.  BISTATIC  SCATTERING 

A  360  degree  pattern  of  the  pressure  amplitude  scattered  from  a  styrofoam 
sphere  was  obtained  to  compare  the  results  with  theoretical  predictions  for  a 
perfectly  reflecting  sphere.  Since  a  sphere  is  perfectly  symmetrical,  it  was  not 
necessary  to  take  the  entire  360  degrees  of  data  but  only  180  degrees  of  it. 

Theoretical  equations  have  been  derived  and  published  by  several  authors  that 
predict  scattering  of  acoustic  plane  waves  from  a  sphere  [Ref.  16,17].  They  tend  to 
be  quite  complicated  functions  involving  sums  of  Hankel  functions.  Figure  41  was 
adapted  from  Medwin  and  Clay  [Ref.  16:pp.  1 92-1 93 J  and  shows  computed 
scattering  patterns  in  relative  amplitude  for  a  perfectly  reflecting  sphere  in  water 
for  several  ka  values.  The  incident  plane  wave  is  arriving  from  the  right.  The 
patterns  show  that  for  increasing  ka  the  major  forward  lobe  gets  larger  in  amplitude 
and  narrower  in  width.  More  structure,  i.e.,  secondary  lobes,  develop  as  ka 
increases. 

Bistatic  data  were  obtained  for  a  ka  =10  using  a  three-inch  diameter  styrofoam 
sphere  and  are  shown  in  Figure  42.  The  incident  wave  arrives  from  the  bottom. 
Frequency  was  61.5  kHz  and  an  eight  cycle  pulse  was  used.  The  graphed  data 
revealed  some  structure  in  the  forward  scattering  region  and  a  large  forward 
scattering  lobe.  Comparison  with  the  theoretical  prediction  for  a  rigid  sphere 
(dashed  line  in  Figure  42)  show  that  the  salient  features  of  the  scattered  pattern 
were  duplicated. 
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Figure  36.  Styrofoam  Sphere  Reflection,  106  kHz 


Figure  40.  Styrofoam  Sphere  Reflection,  106  kHz 
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Figure  42.  Bistatic  Form  Function  of  a  Stvrofoam  Sphere 

ka  =  10 
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A  styrofoam  cylinder  of  diameter  2.5  cm  and  length  20  cm  with  flat  ends  was 
used  for  bistatic  scattering  measurements.  The  particular  diameter  was  chosen 
because  the  styrofoam  was  commercially  available  in  that  size.  The  length  was 
selected  using  two  criteria:  (1)  the  length  was  small  enough  so  that  the  receiver 
could  be  placed  in  the  far  field  of  the  sound  reflected  from  the  cylinder  and  (2)  it 
was  large  enough  to  clearly  present  a  cylindrical  shape. 

Since  the  F41  was  not  omnidirectional,  there  was  concern  of  approximating  a 
plane  wave  at  the  target.  The  lower  the  frequency,  the  better  the  pressure 
distribution  across  the  target  because  of  the  larger  major  lobe,  but  the  wavelength- 
to-target  length  ratio  would  be  less  and  the  signal-to-noise  ratio  also  became  poorer 
at  lower  frequencies.  A  frequency  of  70  kHz  was  chosen  at  a  source-to-target 
distance  of  125  cm.  Figure  43  shows  the  incident  pressure  distribution  as  a  function 
of  distance  off  the  axis  at  a  distance  of  125  cm.  The  lower  curve  is  for  the 
monostatic  case  where  the  receiver  shadowed  some  of  the  incident  energy  with  a 
receiver-to-target  center  distance  of  50  cm.  A  two  percent  pressure  amplitude 
difference  existed  for  the  shadowed  case.  An  eight  percent  decrease  in  pressure 
amplitude  existed  when  not  shadowed. 

The  cylinder  was  mounted  to  a  1/8-inch  diameter  metal  tube  which  had  a  2  X 
1/2  X  1/1 6-inch  flat  plate  epoxied  to  it  to  provide  a  means  of  rotating  the  cylinder. 
The  flat  plate  was  forced  into  the  interior  of  the  cylinder  and  provided  a  snug  fit. 
The  effect  of  the  plate  on  the  reflections  was  assumed  to  be  neglible  since  it  was 
completely  surrounded  by  the  styrofoam.  The  tube  was  as  small  as  possible  to 
minimize  reflections.  A  string  with  weights  was  attached  to  keep  the  styrofoam 
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Figure  43.  Pressure  Distribution  at  70  Khz 
cylinder  submerged.  The  thin  tube  was  soldered  into  a  3/8-nch  diameter  tube 
which  supported  the  entire  arrangement.  The  cylinder  was  level  before  entering  the 
water  and  appeared  to  remain  that  way  when  in  the  water.  Figure  44  illustrates  the 
support  setup . 

Data  were  recorded  for  nine  bistatic  angles.  Time  discrimination  of  the 
reflected  pulse  was  used  whenever  possible.  The  souarce-to-receiver  distance  was 
125  ±  2  cm  in  all  cases  with  the  F41  being  driven  at  177  V  rms.  The  Ithaco 
preamplifier  setting  was  20X  with  a  3-300  kHz  bandwidth.  Receiver- to-cy Under 
center  distance  varied,  but  the  receiver  was  always  in  the  far  field  of  the  cylinder. 
All  data  were  extrapolated  to  one  meter  from  the  target  center  for  comparison.  A 
20  cycle  pulse  was  used. 

The  cylinder  was  initiaUy  aUgned  broadside  to  the  source  by  eye  to  within  an 
uncertainty  of  about  ±  4  degrees.  Assuming  symmetry,  data  were  taken  for  180 
degrees  of  rotation  and  several  points  were  taken  beyond  that  limit.  Figure  45 
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Figure  44.  Support  Platform  Setup 
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illustrates  the  angles  measured  in  the  experiments.  The  varied  angle,  A,  is  the  angle 
made  between  the  cylinder's  normal  and  the  line  connecting  the  center  of  the  target 
to  the  source  axis.  Angle  B  is  the  bistatic  angle. 

Source 
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Figure  45.  Measured  Angles  in  Experiment 
The  results  of  the  experiment  are  shown  in  Figures  46  to  55.  In  most  cases, 
strong  peaks  were  observed  when  angle  A  was  one-half  the  bistatic  angle  and 
secondary  peaks  with  smaller  amplitude  occurred  90  degrees  from  that.  This 
geometry  seems  to  indicate  specular  reflection  off  the  long  cylinder's  surface  and  its 
flat  ends  respectively.  As  the  bistatic  angle  increases  the  secondary  peak  broadens 
and  eventually  disappears.  Well  into  the  forward  scattering  region,  greater  than 
150  degrees  of  bistatic  angle,  the  data  became  less  predictable  The  distinction  of  the 
primary  and  secondary  peaks  became  less  clear  with  additional  peaks  showing  up. 
At  a  bistatic  angle  of  180  degrees  no  distinctive  peaks  were  present  and  amplitude 
became  very  large  in  comparison  to  other  bistatic  angles. 
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vn.  CONCLUSION  AND  RECOMMENDATIONS 

A  data  acquisition  system  using  the  Tektronix  2430  digital  oscilloscope  was 
implemented  to  extract  scattered  waveforms  in  a  noise-filled  water  environment. 
The  system  was  verified  to  work  and  reflected  waveforms  were  collected. 
Preliminary  work  revealed  that  the  relative  movement  of  the  source  and  receiver 
between  collection  of  waveforms  was  critical  in  obtaining  accurate  data  using  a 
digital  sampling  system  with  a  subtraction  routine.  Data  collection  was  slow  and 
tedious. 

Several  aspects  influencing  scattered  waveforms,  including  target  composition, 
frequency,  and  pulse  length  were  observed  and  recorded.  The  effects  of  these  on 
the  waveforms  of  the  reflected  pulse  produce  difficutlies  in  measuring  the  intensity 
of  pulses  with  waveforms  not  identical  of  the  incident  pulse. 

Styrofoam  targets  were  used  to  circumvent  the  problem  of  waveform 
dependence  on  frequency,  pulse  length,  and  composition.  Bistatic  measurements 
were  taken  using  a  sphere  and  a  cylinder  as  targets.  The  data  from  the  sphere 
compared  well  with  theoretical  predictions  for  a  perfectly  reflecting  sphere. 
Readily  available  data  were  not  found  for  comparison  with  the  results  obtained  for 
the  cylinder.  The  goal  of  matching  wavelength-to-target  length  ratios  observed  in 
oceanic  conditions  was  not  done  in  this  work.  The  USRD  type  E27  transducer 
would  be  an  excellent  source  for  this  application  since  it  has  a  frequency  range  100 
to  500  kHz.  Its  small  size  makes  its  beamwidths  very  large  which  is  beneficial  when 
receiver-target  distance  is  limited  by  tank  dimensions.  The  E27  is  also  very 
lightweight  and  easily  moved.  The  effect  of  the  jitter  would  have  to  be  evaluated  at 
the  higher  frequencies. 


65 


The  bistatic  radar  cross-section  theorem  could  not  be  compared  to  the  data 
collected  in  this  report.  The  wavelengths  used  in  the  experiments  are  not  close  to 
the  "vanishing  limit"  as  stated  in  the  theorem.  One  common  trend  between  the 
bistatic  data  and  the  theorem  was  that  the  data  became  less  predictable  as  the  bistatic 
angle  approached  180  degrees. 

For  future  work,  a  physical  structure  must  be  made  to  facilitate  positioning  of 
the  source,  target  and  receiver  in  three  dimensions  so  that  data  collection  can  be 
completely  automated  and  accuracy  improved.  The  rotation  of  supports  should 
have  a  digital  device  to  allow  control  by  the  computer.  The  design  and  construction 
of  such  a  system  would  not  be  a  trivial  project  and  would  require  a  great  amount  of 
ingenuity  and  experience  with  electronic  devices.  The  payoff  of  having  such  a 
system  would  be  great  to  students  studying  underwater  acoustics  and  in  particular  to 
anyone  studying  scattering  effects.  Once  a  positioning  system  is  built  and 
completely  automated,  data  collection  will  become  simpler  and  more  time  can  be 
devoted  to  analyzing  it. 
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APPENDIX  A 

COMPUTER  LISTING  OF  "TEX" 


PROGRAM  TEX 


DATA  AQUISITION  FOR  BI- STATIC  TARGET  STRENGTHS 
USING  THE  TEKTRONIX  2430  DIGITAL  OSCILLOSCOPE 

THIS  PROGRAM  TRANSFERS  WAVEFORM  DATA  FROM  THE  TEX  2430 
TO  THE  COMPUTER  WHERE  IT  IS  CONVERTED  TO  ACTUAL  VOLTAGE 
VALUES  AND  STORED  ON  A  FLEXIBLE  DISK.  INPUT  FROM  THE 
USER  REQUIRES  FILE  NAME.  DESCRIPTION,  AND  CONTINUATION 
RESPONSES.  PREVIOUS  SETUP  OF  THE  TEX  2430  IS  REQUIRED. 

A  THOROUGH  KNOWLEDGE  OF  THE  GPIB  CODES  IS  NECESSARY.  SEE 
THE  OPERATOR'S  MANUAL  AND  THE  INTERFACING  MANUAL  FOR 
THE  TEX  2430. 


INITIALIZATION 


PRINTER  IS  CRT 
OPTION  BASE  0 

DIM  VdataSHSl  .W$r40I  (NAME  OF  DATA  FILE  AND  DESCRIPTION 
DIM  A$f200J .8*150001 .D$I5J f PRE-WAVEFORM  AND  WAVEFORM  STRINGS 
DIM  V< 1023). Tin( 1023), Volt( 1023)!DATA  ARRAYS 
DIM  X*[ 10] ,Ys*( 10) ,Vs*( 101 ,Ts*I 10) 

DIM  An*( 1 1 


TEKTRONIX  SETUP 


OUTPUT  7045”PATH  OFF" ! SHORTENS  THE  DATA  STRING  SENT  BY  THE  2430 


INPUT  FROM  USER 


INPUT  "NAME  OF  THE  DATA  FILE?" . Vdata* 

INPUT  "DOES  FILE  ALREADY  EXIST?! Y- YES)" .An* 

IF  An»-"Y"  THEN  PURGE  Vdata* 

INPUT  "INPUT  FILE  DESCRIPTION?! 40  CHARACTERS  MAX)".U* 

BEEP 

PRINT  "ALLOWING  DATA  DISK  INSERTION.  PRESS  CONTINUE  TO  RESUME.' 
PAUSE 
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500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 


CREATE  BOAT  VdataS . 1 . 17000 
ASSIGN  eOfile  TO  VdataS 
PRINT 


OBTAIN  WAVEFORM 


PRINT  "OBTAIN  DESIRED  WAVEFORM  ONTO  TEXTRONIX" 
PRINT  "PRESS  CONTINUE  TO  RESUME" 

PAUSE 


!  RECEIVING  PRE-WAVEFORM  DATA  AND  THE  WAVEFORM 

iiMtiMtttmffltiitttwttttttitttttttitwttittttwtMtttHiiitutumitiKHtitninwmMtwiwtitmutiiMttimiiiwtt 

! 

OUTPUT  704; "WFM?  XIN’MTIME  INCREMENT  BETWEEN  POINTS 
ENTER  704 ;X$ 

OUTPUT  704; "WFM?  YMU"! VOLTAGE  DIGITIZING  INCREMENT 
ENTER  704 ;Ys* 

OUTPUT  704; "CHI?  VOL" ! VOLTS/DI V 
ENTER  704;Vs$ 

OUTPUT  704; "HOR?  BSE" IB  TIME  SCALE 
ENTER  704;Ts$ 

! 

( 

OUTPUT  704;"CURV?" 

ENTER  704 -.B* 


SEPARATION  OF  ASCII  STRING  INTO  DATA  POINTS 


N=LEN< B*> 

D$»"" 

lnc“0 

! 

FOR  I«1  TO  N 

IF  B«n.I)-","  THEN 

V(Inc)-VAL(DS) 

Inc'Inc+l 

D$-"" 

GOTO  950 
END  IF 

0$*^D$&B$[I.I1 
NEXT  I 
! 

!  CONVERSION  TO  REAL-TIME  DATA 

I  11  n  1MI 11  triMtlWMUMMMMWW 
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1010  '  X$=TIME  INCREMENT  BETWEEN  POINTS  DIV/50 

I0Z0  I  YS*'VOLTS  BETWEEN  POINTS  DIV/2S 

t030  !  VS*-VOLT5/DIV 

1040  I  TS$=B  TRIGGER  TIME  SCALE 

1050  I 

10G0  Del-tin«VAL<X$) 

1070  DelvoI=VAL< Y5$) 

1080  Vscale°VAL< Vs«> 

1090  Tscale-WAL(Ts$) 

1 100  ! 

1110  OUTPUT  eOf ile:WS.Tscale,Vscale 
11Z0  FOR  1-0  TO  1023 
1130  Tin(I)-I»Deltxi-« 

1140  Voltd  >-V(  I  )*Delvol 
1150  NEXT  I 
1  1E0  ! 

1170  PRINT  "WRITING  TO  FILE" 

1180  OUTPUT  eDfile;Tiin(  •)  .Volt(  *) 

1 190  ASSIGN  eOf ile  TO  * 

1200  BEEP 
1210  I 

1230  ! 

1240  INPUT  "DO  YOU  WANT  ANOTHER  WAVEFORM?  Y=YES",An$ 
1250  IF  An$="Y"  THEN  GOTO  430 
12G0  I 

1270  PRINT  "PROGRAM  TERMINATED" 

1280  END 
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APPENDIX  B 

COMPUTER  LISTING  OF  "PLOT 


10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
1 10 
120 
130 
140 
150 
160 
170 
180 
130 
200 
210 
220 
230 
240 
250 
260 
270 
280 
230 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 


PROGRAM:  PLOT 

THIS  PROGRAM  TAKES  TUO  FORMATTED  BOAT  FILES  AND  SUBTRACTS 
ONE  FROM  THE  OTHER.  IT  WAS  DESIGNED  FOR  USE  USING  UAVEFORMS 
OBTAINED  ON  THE  TEKTRONIX  2430  DIGITAL  OSCOPE.  DATA  IS  TAKEN 
FROM  THE  FILES  IN  THE  FOLLOWING  ORDER. 

U$  -  FILE  DESCRIPTION  STRING 

TSCALE  =  TIME  SCALE  OF  THE  SCOPE  ( SEC/DI V) 

VSCALE  =  VOLTAGE  SCALE  OF  THE  SCOPE  ( VOLT/DI  V) 

T(  0-1023)  •=  TIME  ARRAY  (SEC) 

V( 0-1023)  =  VOLTS  ARRAY  (VOLTS) 

THE  SUBTRACTION  RESULTS  ARE  THEN  ABLED  TO  BE  PLOTTED  AS  HARD 
COPY  ON  THE  THINK  JET  OR  PEN  PLOTTER.  A  STANDARD  OPTION  IS 
INCORPORATED  FOR  EASE  UHEN  USING  THE  PROGRAM  FREQUENTLY  UITH 
THE  SAME  PARAMETERS.  THE  PROGRAM  ALLOWS  FLEXIBILITY  IN 
CHANGING  THE  PLOT'S  PARAMTERS  BEFORE  THE  HARD  COPY. 


OPTION  BASE  0 

DIM  X.da-tal  1023) . Y_da-ta(  1023)  ,Y_axis_units$I  32) 

DIM  Head_lbl$lG0] .Y_axis_lbl$lB0] ,X_axis_lbl$[32) 

DIM  X_axis_unit5$[321 .Clear_crt$I21 , C$121 ,Sub_head_lbl$[ 32 ) 
DIM  X_char$t8] .VdataSE 15) .W$[40] ,Vdata2$[ 15) ,Y_char$[Bl 
DIM  XI ( 1023) ,Y)< 1023) .X2(  1023), Y2(  1023) 
Clear_crt$-CHR$(255)iCHR$(75) 

C$«=CHR$(255)&"K" 


GOSUB  Load_file 

INPUT  "DO  YOU  WISH  TO  USE  THE  STANDARD  PARAMETERS?  Y-YES",Dec2$ 
IF  Dec2$-"Y"  THEN 


GOSUB  Stand _parans 
GOTO  420 
END  IF 

GOSUB  Initlal.parans 
A$-"N" 

GOSUB  Plot_data 

INPUT  "WISH  TO  CHANGE  CURRENT  GRAPHICS  PARAMETERS?  Y»YES" .Dec2$ 
IF  Dec2$=”Y"  THEN  GOSUB  Changc.parans 


460  INPUT  "ftNOTHER  DATA  FILE  PLOTTED  ON  SAHE  HARD  COPY  GRAPH?" . Decb$ 
470  WHILE  Oec6S-"Y" 

480  GOSUB  Load_file 

490  INPUT  "LINE  TYPE<  1 -SOLID. 4-OASH. 5-816  DASH.G-DASH-DOT" .Line_nun 

500  INPUT  "PEN  NUMBER  1  TO  8".Pen_nuii 

510  GOSUB  Sawe  graph 

5Z0  GOTO  460 

530  END  WHILE 

540  INPUT  "DO  YOU  WISH  TO  CHANGE  ANY  MORE  PARAMTERS7  Y-YES".Dec2$ 

550  IF  DecZ$-"Y"  THEN  GOTO  450 

560  INPUT  "DO  YOU  WISH  TO  OUW  GRAPHICS?  Y-YES".A$ 

570  IF  A»-"Y"  THEN  GOSUB  Plot.data 

580  INPUT  "DO  YOU  WANT  TO  PLOT  A  DIFFERENT  GRAPH?  Y-YES".Dec7$ 

590  IF  Dec7S-"Y"  THEN  GOTO  350 

600  GOSUB  End_progran 

610  ! 

630  !  STANDARD  PARAMETERS 

650  ! 

660  Stand_parans :  I 
670  Head_lbl$=U$ 

680  SuQ_head_lblS=VdataS 
690  X_char$“"K" 

700  Y_char$“"DD, DD" 

710  Y_axas_lbl$“"Volts" 

720  X_ax xs_lbl*-‘‘Reiat ive  Time" 

730  Y_ax i9_unit6S»"Volts" 

740  X_ax is_unlts*“"Hicroaeconds" 

750  Left-0 

760  Right-19*T8caie 

770  Bot'tom*-V5caie*4 

780  Top-V9cale»4 

790  X_tic_sp-Tscale/Z 

800  Y_tic_sp-V5cale/5 

8»0  X_maj_ct-Z 

820  Y_maj_c't-5 

830  Hai_tic-2 

840  Xgridtic-Tscale 

850  Ygridt ic-V»cal e 

860  Xgr i dmaj - 1 

870  Ygr i dmaj -  1 

880  MajQrid_32-Z 

890  Y_incmt-2»VBcale 

900  X_incmt“4»Tscale 

910  Line_num-1 

9Z0  RETURN 

930  ) 

950  !  INITIAL  PARAMETERS 

960  l«»«»»»0«0«««0«<«<W0«r«<m01WW»««»»»«IOWIIMWRWWt«W»4««B88»>l»l«WWWW 
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I 


970  I 

980  Init ial_param5 : 

990  X_tic_5o=*0 

1000  Y_tic_sp»0 
1010  X_naj_ct»l 
10Z0  Y_naj_ct«l 
1030  Maj_tic“Z 
1040  Xgricliic«0 
1050  Ygridtic»0 
10G0  Xgridnaj=l 
1070  Ygridnaj“l 
1080  Majgrid_5Z»Z 
1090  INPUT  "GRAPH  HEADING  LABEL. 32  CHAR.  MAX" .Head_lbl$ 

1100  INPUT  "GRAPH  SUB  LABEL. 3Z  CHAR.  MAX" . Sub_head_l bl $ 

1110  INPUT  "IMAGE  FOR  X  AXIS  NUMERICS.  8  CHAR.  MAX" . X_char$ 
1120  INPUT  "IMAGE  FOP  Y  AXIS  NUMERICS.  8  CHAR.  MAX".Y_char$ 
1130  INPUT  "Y  AXIS  L--'HEL.32  CHAR.  MAX"  .  Y_ax  1 5_lbl$ 

:i40  INPUT  "X  AXIS  LABEL, 3Z  CHAR.  MAX" ,X_ax is_lbl * 

M50  INPUT  "Y  AXIS  UNITS. 32  CHAR.  MAX" .Y_axis_units$ 

M60  IN>^UT  "X  AXIS  UNITS.3Z  CHAR.  MAX"  . X_ax i s.unitsS 

1170  DI-P  "UINOOU  LOOKED  AT .LEFT . RIGHT . BOTTOM , TOP" t 

1171  INPUT  Lef t , Riyht , Bottom , Top 

1180  INPUT  "X  AXIS  TICK  SPACING .0EF=0(  NO  TI CKS > " . X_t ic_sp 

1190  INPUT  "Y  AXIS  TICK  SPACIN6.DEF-0(  NO  TICKS )", Y_t ic_sp 

1200  INPUT  "NUM  TIC  INTVLS  BTUN  MAJ  TICS  X-DIR,OEF-l" .X_naj_ct 
1210  INPUT  "NUM  TIC  INTVLS  BTUN  MAJ  TICS  Y-DIR,DEF-1" .Y_«aj_ct 
1220  INPUT  "SIZE  OF  MAJOR  TICKS. DEF-2" .Maj_tic 
1230  INPUT  "GRID  X-TIC  SPACING ,OEF-0< NO  TICKS)" .Xgridtic 
1240  INPUT  "GRID  Y-TIC  SPACING.DEF-0( NO  TICKS)" .Ygridtic 
1250  INPUT  "NUM  TICKS  BTUN  MAJOR  TICKS  Y-DIR, DEF» ) " . Ygr idnaj 
1ZG0  INPUT  "NUM  TICKS  BTUN  MAJOR  TICKS  X-DIR.DEF-1"  .Xgndmaj 
1270  INPUT  "SIZE  OF  MAJOR  GRID  TICKS. DEF-2" .MaJgrid.sz 
1280  INPUT  "Y-AXIS  INCREMENT"  .Y_inciit 

1290  INPUT  "X-AXIS  INCREMENT" ,X_xnc«t 

1300  DISP  "LINE  TYPE( 1 -SOLID. 4-DASH. 5-BIG  DASH.G-OASH-OOT)"j 

1301  INPUT  Line_nun 

1310  > 

1 320  RETURN 

1 330  !  _  _ 

1350  1  LOAD  DATA  FILE 

1  360  I  «l1«1t<1ltW«g»IW4Mt<mi«W1<»W1l1MtllM«»«tt«lt«»«1>«1)<»»MMMI1t11»«<tt««ir 
1370  ! 

1380  Load_file:  ! 

1 390  ! 

1400  INPUT  "INPUT  NAME  OF  ORTA  TILE  r'.Vdata* 

1410  INPUT  "INPUT  NAME  OF  DATA  FILE  2  ".Vdata2* 

1420  ASSIGN  eOfI  TO  Vdata* 
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1430 
1440 
14^0 
1460 
1470 
1480 
1490 
1S00 
1510 
I5Z0 
1530 
1540 
1550 
1560 
1570 

1580  Plot_data;  ! 

1590  ! 

1600  IF  8$-"Y"  THEN 

16)0  OUTPUT  KBO;Clear_crt*; 

I6Z0  GINIT 

1630  PLOTTER  IS  CRT , “INTERNAL" 

1640  GOTO  1770 

1650  END  IF 

1660  INPUT  "IS  THIS  TO  BE  ft  HARDCOPY  PLOT?  Y-YES",Decl$ 
1670  IF  D«‘cl$-"Y"  THEN 
1680  GfNir 

1690  PLOTTER  IS  705."HPGL" 

17<>^  GOTO  1770 

1710  ELSE 

17:0  OUTPUT  KBO;Clear_cr-tSj 

1730  GINIT 

1740  PLOTTER  IS  CRT INTERNAL" 

1750  GOTO  1770 

1760  END  IF 
1770  GRAPHICS  ON 
1780  PEN  1 

1790  X_gdu_«ax-100*MAX( 1 .RATIO) 

1800  Y_Qdu_«ax-100*nAX(1 ,1 /RATIO) 

1810  LORG  5 
1820  CSIZE  4 

1830  MOVE  .48* K  odu  i»iax-.93*Y  adu  wax 
1840  LABEL  Head_lbl* 

1850  MOVE  . 48*X  ndu  wax..89*Y  odu  wax 

1860  CSIZE  3 

1870  LABEL  Sub_head_lbl$ 

1 880  DEG 
1890  LOIR  90 
1900  CSIZE  4 


ENTER  eOf 1 ? W$ . Tscal e . Vscal e,X1(*>.Yl(*) 
ASSIGN  eOfl  TO  • 

ASSIGN  @Df2  TO  VdataZ$ 

ENTER  eOf ZjUS. Tscal e. Vscal e.X2(  * ) , YZ(  » ) 
ASSIGN  eOfZ  TO  • 

MAT  X_data-  XI 
MAT  Y_data«  Y1-Y2 
PRINT  "DATA  FILES  READ" 

BEEP 

RETURN 

I 


PLOT  DATA 
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1910 
1920 
1930 
1940 
1950 
1950 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
21  10 
2  120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2710 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2  330 
2340 
2  350 
2350 
2  370 
2  380 


MOVE  . 04»X  qdii  nax..55»Y  Qdu  wax 
LABEL  Y_axi5_lbl$ 

LORG  5 
LDIR  0 

MOVE  . 48*X_gdu_max . . I5»Y  qdu  wax 
C5IZE  4 

LABEL  X_axis_lbl$ 

MOVE  . 48» X_Qdu_max , . I  9* Y_gdu_nax 
CSIZE  3 

LABEL  ;  X_ax  is_Linits#; 

MOVE  . 70*X  qdu  wax..74*Y  qdu  wax 

VIEWPORT  .  15»X_gdu_Max  , .  85»X_Qdu_rtax  , .  25»Y_gdu_inax  , .  86»Y_Qdu_aax 
UINOOU  Lef t .Right .Bot ton, Top 

AXES  X_t  ic_sp ,  Y_t  lc_5p  ,Lef  t , Bottom ,  X_naj_ct ,  Y_piaj_ct , Maj _t  i c 
AXES  X_t ic_sp , Y_t ic_sp .Right . Top . X_maj _ct , Y_maj _ct ,Maj_t ic 
•  GRID  Xgr idt  ic  .  Ygr idt  1C ,Lef  t  .Bottom .  XgndmaJ  .Ygridmaj  .Majgrid_5z 
CLIP  OFF 
CSIZE  3.-75 
LORG  6 

FOR  I=Left  TO  Right  STEP  X_incmt ! LABELING  OF  X  AXIS 
MOVE  I .( Bottom- .04»( Top) > 

LABEL  USING  X_charS ; I / 1 . E-G 
NEXT  1 
LORG  8 

FOR  I»Bottom  TO  Top  STEP  Y_incmt  •  LABELI NG  OF  Y  AXIS- 
MOVE  (Left-,4»TscaU)  .1 
LABEL  USING  Y_char$5l 
NEXT  I 
PENUP 
PEN  2 

LINE  TYPE  Line_num 
CLIP  ON 

FOR  1-0  TO  t023'PLOTTING  OF  DATA 
PLOT  X_data< I ) . Y_data( I ) 

NEXT  I 
PENUP 
CLIP  OFF 

IF  A$="Y"  THEN  DUMP  GRAPHICS  #701 
RETURN 


Same  oraoh'  ! 

PEN  Pen_num 
LINE  TYPE  Line_num 
CLIP  ON 
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Z390  FOR  1=0  ro  I0i;3 

Z400  PLOT  X_data( I ) . Y_data( I  ) 

Z4t0  NEXT  I 
Z4Z0  CLIP  OFF 
Z430  PENUP 
Z440  RETURN 
2450  I 

Z470  I  CHANGE  GRAPHICS  PARAMETERS 

2490  I 

2500  Change _para«5:  > 

2510  OUTPUT  XBD:Cloar_crt*; 

2520  PRINT  "SELECT  FROM  THE  FOLLOWING  TO  CHANGE" 

2530  PRINT  "CASE  1:  GRAPH  HEADING" 

2540  PRINT  "CASE  2:  Y  AXIS  LABEL" 

Z550  PRINT  "CASE  3:  X  AXIS  LABEL" 

2SS0  PRINT  "CASE  4:  Y  AXIS  UNITS" 

2570  PRINT  "CASE  5:  X  AXIS  UNITS" 

2580  PRINT  "CASE  G:  WINDOW  UNITS" 

2590  PRINT  "CASE  7:  X  AXIS  TICK  SPACING" 

2G00  PRINT  "CASE  Q:  Y  AXIS  TICK  SPACING" 

2G10  PRINT  "CASE  9:  TICK  INTVL  BTWN  MAJOR  TICKS  IN  X  DIRECTION" 

ZG20  PRINT  "CASE  10:  TICK  INTVL  BTWN  MAJOR  TICKS  IN  Y  DIRECTION" 

2B30  PRINT  "CASE  11:  SIZE  OF  MAJOR  TICKS" 

ZB40  PRINT  "CASE  12:  X  TICK  GRID  SPACING" 

2G50  PRINT  "CASE  13:  Y  TICK  GRID  SPACING" 

ZGG0  PRINT  "CASE  14:  GRID  TICK  INTVL  BTWN  MAJOR  TICKS  IN  X  DIRECTION" 

2B70  PCINT  "CASE  15:  GRID  TICK  INTVL  BTWN  MAJOR  TICKS  IN  Y  DIRECTION" 

2680  PRINT  "CASE  IB:  MAJOR  GRID  TICK  SIZE" 

ZB90  PRINT  "CASE  17:  Y  INCREMENT" 

2700  PRINT  "CASE  18:  X  INCREMENT" 

2710  PRINT  "CASE  19:  LINE  TYPE" 

27Z0  PRINT  "CASE  20:  GRAPH  HEADING  SUB  LABEL" 

2730  PRINT  "CASE  Zl:  X  AXIS  NUMERIC  IMAGE" 

2740  PRINT  "CASE  22:  AXIS  NUMERIC  IMAGE" 

2750  1 

2760  INPUT  "DO  YOU  DESIRE  TO  CHANGE  ANY  PARAMETERS?  Y-YES".Dec4* 

2770  IF  Dec4S-"Y"  THEN 

2780  INPUT  "ENTER  PARAMETER  NUMBER" .Change.nuB 

2790  ELSE 

2800  GOT)  35G0 

2810  END  IF 

2820  SELECT  Change_nu« 

2830  CASE  1 

2840  PRINT  "CURRENT  HEADING  LABEL" .Head_lbl» 

2850  INPUT  "NEW  HEADING  LABEL" .Head_lbl$ 

28G0  CASE  Z 
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Z870  PRINT  "CURRENT  Y  AXIS  LABEL" . Y_ax i5_lbl$ 

Z880  INPUT  "NEU  Y  AXIS  LABEL" . Y_ax i 5_1 bl $ 

Z890  CASE  3 

Z900  PRINT  "CURRENT  X  AXIS  LABEL" , X_ax is_l bl * 

Z310  IW>UT  "NEW  X  AXIS  LABEL" . X_ax is_l bl S 

Z9Z0  CASE  4 

Z930  PRINT  "CURRENT  Y  AXIS  UNI TS" . Y_ax is_uni tsS 

Z940  INPUT  "NEU  Y  AXIS  UNITS" .Y_axis_unit5* 

Z950  CASE  5 

Z9G0  PRINT  "CURRENT  X  AXIS  UNITS" .X_axis_unit5t 

Z970  INPUT  "NEU  X  AXIS  UNITS" .X_ax i6_units$ 

Z380  CASE  G 

Z990  PRINT  "CURRENT  UINDOU" .Left .Right .Botion . Top 

3000  INPUT  "NEU  UINDOU" .Left .Right .Bottom, Top 

3010  CASE  7 

30Z0  PRINT  "CURRENT  X  AXIS  TIC  SPACING" . X_t ic_sp 

3030  INPUT  "NEU  X  AXIS  TIC  SPACING" .X_tic_sp 

3040  CASE  8 

3050  PRINT  "CURRENT  Y  AXIS  TICK  SPACING" .Y_tic_5p 

3060  INPUT  "NEU  Y  AXIS  TIC  SPACING" . Y_t ic_sp 

3070  CASE  9 

3080  PRINT  "CURRENT  TIC  INTVL  BTUN  MAJ  TICS  X  01 R" . X_maj _ct 

3090  INPUT  "NEU  TICK  INTVL  BTUN  hAJOR  TICKS  X  DIR"  .X_ntaj_ct 

3100  CASE  10 

3110  PRINT  "CURRENT  TICK  INTVL  BTUN  MAJOR  TICKS  Y  DIR" .Y_maJ_ct 

31Z0  INPUT  "NEU  TICK  INTVL  BTUN  MAJOR  TICKS  Y  DIR" . Y_ma]_ct 

3130  CASE  n 

3140  PRINT  "CURRENT  MAJOR  TICK  SIZE" .Maj_tic 

3150  INPUT  "NEU  MAJOR  TICK  SIZE" .Maj_tic 

31B0  CASE  IZ 

3170  PRINT  "CURRENT  X  TICK  GRID  SPACING" ,Xjgrid_tic 

3180  INPUT  "NEU  X  TICK  GRID  SPACING" ,X_orid_tic 

3130  CASE  13 

3Z00  PRINT  "CURRENT  Y  TICK  GRID  SPACING" ,Y_grid_tic 

3Z10  INPUT  "NEU  Y  TICK  GRID  SPACING" .Y_grid_tlc 

32Z0  CASE  14 

3230  PRINT  "CURRENT  GRID  TIC  INTVL  BTUN  MAJ  TICS  X  DIR" .X_Qrid_naj 

3240  INPUT  "NEU  GRID  TIC  INTVL  BTUN  MAJ  TIC  X  DIR"  .X_orid_i«iaj 

3250  CASE  15 

32G0  PRINT  "CURRENT  GRID  TIC  INTVL  BTUN  MAJ  TICS  Y  DIR" . Y^jjrld.naj 

3270  INPUT  "NEU  GRID  TIC  INTVL  BTUN  MAJ  TIC  Y  DIR" . Y_Qrid_naj 

3280  CASE  IG 

3290  PRINT  "CURRENT  MAJOR  TICK  SIZE" .Majjgrid.si 

3300  INPUT  "NEU  MAJOR  TICK  SIZE"  .Maj_or-id_5Z 

3310  CASE  17 

33Z0  PRINT  "CURRENT  Y  INCREMENT" . Y_incmt 

3330  INPUT  "NEU  Y  INCREMENT" .Y_incmt 

3340  CASE  18 
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3350  PRINT  ‘CURRENT  X  INCREMENT" . X_incnt 

3360  PRINT  "NEU  X  INCREMENT" .X.incnt 

3370  CASE  19 

3380  PRINT  "CURRENT  LINE  TYPE" .Line_nun 

3390  DISP  "NEU  LINE  TYPE!  1  »S0LI0.4-0flSH  .5-BIG  DASH.G»DftSH-DOT'‘ j 

3391  INPUT  Line_rujn 

3400  CASE  Z0 

3410  PRINT  "CURRENT  HEADING  SUB  LABEL" ,Sub_head_lbl$ 

34Z0  INPUT  "NEU  HEADING  SUB  LABEL" . Sub_head_l bl $ 

3430  CASE  Z1 

3440  PRINT  "CURRENT  X  AXIS  NUMERIC  IMAGE" . X_char* 

3450  INPUT  "NEU  X  AXIS  NUMERIC  IMAGE" . X_charS 

3460  CASE  ZZ 

3470  PRINT  "CURRENT  Y  AXIS  NUMERIC  IMAGE" .Y_char* 

3480  INPUT  "NEW  Y  AXIS  NUMERIC  IMAGE" . Y_ char* 

3490  END  SELECT 

3500  INPUT  "DO  YOU  UISH  TO  CHANGE  ANOTHER  PARAMETER?  Y-YES" .0ec5* 

3510  IF  Dec5*-"Y"  THEN 
35Z0  GOTO  Z510 

3530  ELSE 

3540  GOSUB  PIot_data 

3550  END  IF 
3560  RETURN 
3570  I 

3590  !  ENO  ThC  PROGRAM 

3610  I 

36Z0  End _proQrai<i:  I 

3630  GRAPHICS  OFF 

3640  PRINT  "END  OF  PROGRAM" 

3650  END 
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